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Dust storms include particulate matter that is transported over land and sea with biota that could impact 
downwind ecosystems. In addition to the physico-chemical compositions, organismal diversities of dust 
from two storm events in southern Israel, December 2012 (Evl2) and January 2013 (Evl3), were determined 
by pyro-sequencing using primers universal to 16S and 18S rRNA genes and compared. The 
bio-assemblages in the collected dust samples were affiliated with scores of different taxa. Distinct patterns 
of richness and diversity of the two events were influenced by the origins of the air masses: Evl3 was rich 
with reads affiliated to Betaproteobacteria and Embryophyta, consistent with a European origin. Evl2, 
originated in north-Africa, contained significandy more of the Actinobacteria and fungi, without conifers. 
The abundance of bacterial and eukaryotic reads demonstrates dissemination of biological material in dust 
that may impose health hazards of pathogens and allergens, and influence vegetation migration throughout 
the world. 
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Dust storms, considered major contributors to global aerosols', transport desert soUs to the atmosphere and 
substantially impact the global environment. Estimates of global dust emissions from soUs to the atmo- 
sphere vary between 1 and 3 pg per year^. Climatic variations affect the land surfaces of aeolian (wind- 
driven) systems^, and droughts have significantly contributed to increased dust emissions'''\ The contribution of 
anthropogenic soils, e.g. agricultural fields and industrial yards, are of increasing concern*". Up to 50% of the total 
atmospheric dust originates from disturbed soUs^, which may contain different particle compositions from that of 
natural dust. The windblown dust can travel tens of thousands of kilometers before being deposited', depending 
on the particle characteristics (size, chemistry) and the air-mass properties (e.g., velocity, density, height)". The 
size and chemical compositions vary in space and time and may determine the potential impact on air quality and 
human health' '". 

Airborne dust particles contain different microorganisms that are exceedingly mobile in space and time" '^. For 
example, the number of culturable microorganisms in the US Virgin Islands increased during African dust-storm 
events by about 8 fold, from 0.013 L~' of air under clear atmospheric conditions to O.IOS"''"". Dust storms 
originating in the Saharan desert heavily affect the south-eastern Mediterranean basin, mostly during the winter 
and spring'^; qualitative and quantitative analyses of dust-associated fungal communities revealed distinct pat- 
tern of distribution in the atmosphere of Haifa (Israel)"" '^. 

Airborne dust contains a variety of chemicals and microbial agents such as bacteria, fungi, and viruses where 
some of them are pathogenic and pose a risk to the ecosystem and human health as the clouds traverse 
regions" '^ '". The composition of microorganisms is stiU not well-defined, and taxonomic studies of organisms' 
diversity in the outdoor air have just started to emerge. In recent years, conventional molecular approaches have 
widely been used to study the diversity and community composition of prokaryotes and eukaryotes in air and 
outdoor dusts'^ ''^ but the methods exploited are incomplete. Clone libraries. Restriction Fragment Length 
Polymorphism, Denaturative Gradient Gel Electrophoresis, Ribosomal Intergenic Spacer Analysis and PCR- 
single Strand Conformation Polymorphism are usually restricted to less than 500 sequences or patterns. These 
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approaches, albeit providing general information on the structure of 
explored communities, are therefore not sufficient for meaningful 
comparisons. 

In this study, the high throughput sequencing method (pyrose- 
quencing) was used to explore, the structure of the prokaryotic and 
eukaryotic communities in airborne samples of dust collected in an 
urban environment in the Negev (southern Israel) following two 
distinct desert storm events. 

Results and Discussion 

The arena. Samples collected in an urban environment in the Negev 
(southern Israel) were used as a case study. The Negev is located at 
the margin of the natural dust sources and is frequently subjected to 
such storms with common duration of several hours to one-day^^'^'. 
Recent studies relied mainly on satellite images (MODIS) and the 
HYSPLIT model of air mass backward trajectories to indicate the 
geographic source of the observed bio samples''^' Satellite images 
of the studied region during the dust storm events (each lasted less 
than a day), the air mass trajectories and wind directions at the level 
of the measurement point (consistent with the HYSPLIT trajectories) 
are displayed in Fig. 1. Both storms were associated with a synoptic 
system of the (Cyprus) cold low-pressure that typically moves 
eastward over the Mediterranean Sea^'. The trajectory of the 
December air mass (Evl2) extended over North Africa. Due to a 
deeper cold low in January 2013 (996 mb) than in December 2012 
(1004 mb) the origin of the January air mass (Evl3) was attributed to 
South Europe (Fig. 1). 

Physical and chemical properties. As is generally accepted, stronger 
winds and lower temperatures associated with lower pressure were 
observed here too (Fig. 1): the levels of major air gas pollutants were 



not affected by the dust events and remained relatively low as in non- 
dusty days in the studied region, but those of PMio (particulate 
matter £ 10 [im in diameter) in both storms were about 20-fold 
higher than the background value in the area at non-dusty days 
(42 |ig m"''), as also observed in other strong dust storms in the 
Negev during the last decade'". High levels of PM during dust 
storm events raise air pollution above the standard values of air 
quality'"'^'. In arid environments, hourly PM concentrations 
during dust storms can reach even 10 mg m""* 

The particle size distributions of the dust were typically bi-modal 
for both samples with a cutoff at about 25 |im (Fig. 2). The January 
sample (Evl3) included higher content of finer particles (45.1%) than 
that (34.0%) of December's (Evl2) and smaller sizes of the coarser 
population (peaks at 50 and 60 |im respectively). The possible 
sources of the finer particles, which can be carried longer distances, 
are North- African and South-European soils (see mass trajectories 
and compatible wind directions 10 m above the surface in Fig. 1), 
whereas the possible proxy source is the Negev soils located close to 
Be'er Sheva. In both Evl2 and Evl3 samples, there were relatively 
high (> 22%) frequencies of the class weight of respirable particles 
(< 10 nm) (Fig. 2) that affect human health'. 

Typical elements were found in both dust samples (Table SI), the 
most common of which were Si, as in mineral soils, silt and sand 
fractions (> 50 \imy^ and Ca originating in sedimentary environ- 
ments (such as dry lakes), calcareous soils and rocks in the arid land 
of the Mediterranean basin. The higher percentage of Si in Evl2 
(Table SI) may indicate a longer terrestrial transport of the dust 
(Fig. 1). No significant differences in the contents of Al, Fe, K, Mg 
and Na were found between the dust samples, but some differences 
were found in the contents of minor elements: Co was lacking in 
Evl2; Cu, Os, Ni and Ga were not observed in Evl3. 
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Figure 1 | Satellite images (MODIS) of the studied region during the dust storms along with air mass transport at different heights above ground level 
(AGL) derived from Backward Trajectories model NOAA/ARL HYSPLIT-4 (credit to: www.noaa.gov) and equivalent wind directions at 10 m above 
surface level from the HCMR POSEIDON System (credit to: www.poseidon.hcmr.gr). Red arrowheads indicate the sampling site. Right hand side 
panels: daily recorded averages of major meteorological variables and pollutants. AT - air temperature; RH - relative humidity; WS - wind speed; 
PM - particulate matter; ST - settled dust. 
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Figure 2 | Particle size distribution of the dust samples by a high-resolution laser diffractometer technique (ANALYSETTE 22 MicroTec Plus), with 
the statistical parameters. 



Richness and diversity of organisms in storm dust. The four 
libraries represent true diversity reasonably well: the average 
Good's coverage of the reads in both samples (at 97% cut-off) was 
94.5% for eukaryotes and 86% for bacteria (Table 1). On the other 
hand, the relatively low number of observed Operational Taxonomic 
Units (OTUs; Table 1 and Fig. SI), reaching 46% and 58% saturation 
(Sobs/Schao; Table 1) for eukaryotes and bacteria respectively, 
suggests an under-sampling of the dust; diversity may however be 
biased by tendency of Chao index to overestimate species richness'''. 
The bacterial diversity in both samples was significantly higher than 
the eukaryotic, and the total richness of the dust sample Evl2 was 
higher than of Evl3 (Table 1). 

Venn diagrams (Fig. S2) demonstrate a higher number of bacterial 
and eukaryotic OTUs that were unique in both dust samples, imply- 
ing different origins of the biota. The bacterial samples shared 200 
OTUs (Fig. S2A), representing the majority of the reads, 54% and 
62% for Evl2Bac and Evl3Bac respectively. Most of the unique 1,014 
OTUs of Evl2Bac (Table S2) are singletons (i.e., OTU containing a 
single sequence), implying a relatively low coverage (83%) and high 
values of richness and diversity of this sample. PGR bias associated 
with initial amplicon generation may impose distortions in the 
observed community structure, particularly over-estimating rare 
taxa. Detection limit of the latter may be affected by the approach 
employed in quality control as well, but pyrosequencing itself appears 
not to impose significant bias of overall community structure esti- 
mates'^. Eukaryotic samples shared only 34 OTUs (Fig. S2B, Tables 
S3 and S4), which also represented the majority of the sequences, 
66% and 47% for Evl2Euk and Evl3Euk. Higher diversity and rich- 
ness was also observed in Evl2Euk than in Evl3Euk (Table 1). 

Prokaryotic sequences (Fig. 3A). Reads of 16S rRNA genes that 
were retrieved from both dust storm samples belong to the 



following phyla: Acidobacteria, Actinobacteria, Bacteroidetes, 
Deinococcus-Thermus, Firmicutes, Proteobacteria, Cyanobacteria, 
in addition to unclassified bacteria. Reads of Gemmatimonadetes 
and Chloroflexi were found only in Evl2Bac, and Verrucomicrobia, 
only in Evl3Bac. On a class level, the distributions were similar, with 
the most dominant Alphaproteobacteria (44.0% and 46.7% for 
Evl2Bac and Evl3Bac respectively), Actinobacteria (23.9% and 
17.5%) and Betaproteobacteria (4.0% and 15.5%). 

Actinobacteria. This one of the dominant bacterial phylum/class is 
made up of gram-positive bacteria that are widely distributed in both 
terrestrial and aquatic ecosystems. Some of the Actinobacteria are 
spore forming which range from motile zoospores to specialized 
propagules which resist desiccation and mild heat. Most of the 
sequences are affiliated within the order Actinomycetales (95.7% 
and 99.3% for Evl2Bac and Evl3Bac respectively). The 
Geodermatophilaceae (22.9% and 15.2%), Micrococcaceae (22.9% 
and 42.3%) and unclassified-actinomycetales (27.6% and 23.5%) 
are the dominated families. The Actinobacteria reads belong to 3 
major genera, Kocuria, Arthrobacter and Blastococcus, the species 
of which are found in various environments such as plants, soils, 
sediments and rhizomes. Among all bacterial reads, 4.7% in 
Evl2Bac and 20.5% in Evl3Bac displayed variable similarity (84- 
100%) to the Kocuria species, some with 99% homology to K. rosea 
(KG689298, FJ745378). Several groups are using metagenomic 
framework to study airborne microbial communities""'^. During a 
severe smog event over Beijing, for example, Cao et al.^'' found that 
the most abundant bacterial phylum, order and species were respect- 
ively Actinobacteria, Actinomycetales and Geodermatophilus 
obscurus, Modestobacter marinus, Blastococcus saxobsidens, 
Kocuria rhizophila, and Micrococcus luteus. 



Table 1 | Alpha-diversity indices (97%) based on 454-pyrosequencing data from the dust samples: coverage, sobs (# of OTUs), invSimpson 



and chao parameters" 
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Figure 3 | Distributions of rRNA genes reads, retrieved from the dust samples Evl2 and Evl3. (a) Bacterial 16S reads at the class level. Affiliation of 
classes to different phyla are: 1-4 Acidobacteria; 5 - Actinobacteria; 6-7 - Armatimonadetes; 8-10 - Bacteroidetes; 11 - Chloroflexi; 12 - Deinococcus- 
Thermus; 13-15 - Firmecutes; 16 - Gemmatimonadetes; 17-21 - Proteobacteria; 22 - Verrucomicrobia; 23 - unclassified, (b) Eukaryotic 18S reads at the 
kingdom-class levels. Affiliation of taxonomic categories to highest taxonomic ranks are: 1-16 - Fungi kingdom ( 1-7 - Ascomycota phylum; 8-14 - 
Basidiomycota phylum; 15 - unclassified Dikarya superphylum; 16 - Basal fungal lineages phylum); 17-19 - Metazoa kingdom (17-18 - Arthropoda 
phylum; 19 - unclassified Metazoa); 20-21 - Alveolata superphylum (Ciliophora phylum); 22-28 - Viridiplantae phylum; 29 - unclassified eukaryota 
phylum. 



Kocuria spp. belonging to Micrococcaceae are considered as non- 
pathogenic commensals that colonize the oropharynx, skin and 
mucosa, but K. rosea is an opportunistic pathogen in immune-com- 
promised patients'"' and K. kristinae is associated with acute chole- 
cystitis''. In Evl2Bac and EvlSBac, 11.4% and 14% respectively are 
closely related to Arthrobacter spp. (e.g., HF585203, JX164047, 
JN680244, AJ576068, AB522428), which are widely distributed in 
nature, particularly soU. Their overall pathogenic potential is rather 
low, though some Arthrobacter stains were isolated from clinical 
specimens'"'. Most of the 10.9% and 7.6% of the reads from both 
samples respectively, affiliated with Blastococcus spp., are identical 
to B. saxobsidens (EU977834) and B. aggregatus (FR865889), found 
on rock surfaces in the Mediterranean basin'". Very few sequences 
(11 and 7 respectively) were identical to those of the genus 
Microbacterium found in soil and sewage, some of which have been 
recognized as pathogens in humans, causing e.g., endophthalmitis'*^. 

Alphaproteobacteria. This class is a group of gram-negative bacteria 
that comprise most phototrophic genera, symbionts of plants and 
animals and a group of pathogens. Orders Rhodobacterales (28.9% 
and 59.5%) mdRhizobiales (44.6% and 17.1%), are among the largest 
revealed in both Evl2Bac and Evl3Bac respectively; the former are 
dominant and ubiquitous primary surface colonizers in temperate 
coastal waters of the world"" and the latter include nitrogen fixing 
symbionts and pathogenic to animals and plants'". Alphaproteo- 
bacteria in both samples showed similar distribution on a family 
and genus levels and was dominated by the following genera: 
Paracoccus (5.3% and 15.8%), RubelUmicrobium (15.5% and 
24.5%) and Skermanella (6.7% and 9.3%). Some of the ubiquitous 
spp. of the genus Brevundimonas are considered to be opportunistic 
pathogens in immune-compromised hosts'*^, e.g., B. vesicularis and 
B. diminuta are frequently isolated species in human infections. Only 
16 and 14 reads here were identical to B. vesicularis (KC494336, 
GU430201, HM755555). 

Betaproteobacteria. This is the 3rd most dominant class, with 89% 
and 83% of the reads for Evl2Bac and Evl3Bac respectively, belong- 
ing to the order Burkholderiales that includes several pathogenic 



genera e.g., Burkholderia and Bordetella, but none of these reads 
belongs to them. Among the reads of this order, 43.6% and 62.7% 
are identical to Massilia yuzhufengensis (JQ409016) isolated from ice 
drilled in Yuzhufeng Glacier, Tibetan Plateau, China'"'. 

Unclassified bacteria. A large portion of these reads in Evl2Bac and 
Evl3Bac (19.1% and 8.6% respectively) were revealed by the RDP 
database; when analyzed against NCBI GenBank, 56.5% and 11.9% 
of these were identical to the cyanobacterial species Oscillatoria 
nigroviridis PCC 7112 (CP003614). 

Viability of pathogenic microorganisms. Viability is a meaningful 
parameter in developing predictive models for disease dispersaP'. 
An average of about 70% viable bacteria has recently been estimated 
by fluorescent staining in anticyclone air""*, but the majority (>99%) 
of microorganisms from the environment resist cultivation in the 
laboratory"" ''^, hence culture dependent methods cannot be used to 
assess viability of the pathogens in dust samples. 

Eukaryotic sequences (Fig. 3B). Dothideomycetes. This is the largest 
and most diverse class of ascomycetes that includes several plant 
pathogens (e.g., Phaeosphaeria nodorum, Venturia inaequalis). 
Evl2Euk and Evl3Euk contained 13.5% and 3.15% such reads 
respectively, most of them (8.9% and 1.5%) are highly similar (99- 
100%) to Cladosporium spp. (e.g. JX273066, JN974018, JN546118), 
some of which cause infections of the skin and toenails, sinusitis and 
pulmonary infections. Their airborne spores are significant aller- 
gens'"' and in large amounts they can severely affect asthmatics'". 
Dominant sequences retrieved from indoor dust of urban area in 
central Finland were identical to C. cladosporioides and C. herba- 
rum^^. The rest of the Dothideomycetes reads (4.3% and 1.3%) were 
identical to Alternaria spp. (e.g., KC584600, KC584596, KC584599), 
major plant pathogens and common humans allergens'"*, causing hay 
fever or hypersensitivity reactions that sometimes lead to asthma'". 
Many health disorders are caused by these fungi, which grow on skin, 
mucous membranes, the eyeballs and the respiratory tract^". Alter- 
naria and Epicoccum spp. (producing multicellular dictyospores > 
10 |.im) are abundant allergenic fungi, mostly in the larger particle- 
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size ranges (d^ > 44.7 [im), whereas Cladosporium is an abundant 
allergenic fungal genus, distributed evenly across all the particle-sizes 
ranges'''. Viable microbial populations, including presumptive plant 
pathogens Alternaria infectoria and Chaetomium globosum, have 
recently been detected in Asian air samples even after traveling 10 
days across the Pacific Ocean in the free troposphere, information 
that has significant implications for epidemiology^'^. 

Tremellomycetes. This fugal class {Agaricomycotina, Basidiomycota) 
is a dimorphic, nutritionally heterogeneous group comprising of 
saprotrophs, animal parasites, severe human pathogens and fungi- 
colous species^^. Evl2Euk and Evl3Euk contain 12% and 3.3% of 
such reads that are highly similar to Cryptococcus spp. and 
Filobasidium spp.; about 60% of them are identical to C. albidus 
(HQ231895) that occasionally causes moderate-to-severe diseases, 
specifically meningitis, in patients with compromised immunity^"". 

Mitosporic Ascomycota. These fungi comprise a heterogeneous 
group and represent more than half of Ascomycota lacking a sexual 
state; many pathogenic fungi in plants and mammals, including 
humans, belong to this group""*. Evl2Euk and Evl3Euk contain 
3.9% and 0.5% reads respectively belonging to this group. 

Embryophyta. This is the most familiar subkingdom of green plants 
(Viridiplantae), informally called land plants, excluding green 
algae""^, and was represented by 10.4% and 79.2% in Evl2Euk and 
Evl3Euk respectively. This high percent of plant reads in the latter is 
consistent with southern Europe as the origin of the air mass in 
January 2013. Among them, 3.6% and 2.7% in Evl2Euk and 
Evl3Euk respectively displayed 98-100% homology to leafy trees 
from the Morus spp. (GU476477 and L24398), Moringa oleifera 
(U42786), Plocosperma buxifolium (HQ384684), Metteniusa tess- 
manniana (AM421127) and Oka europaea (L49289). Additional 
3% of Embryophyta reads in each Evl2Euk and Evl 3Euk are identical 
to Bryum pseudotriquetrum (KC291525), Blindia acuta (AF023681) 
and Pottia truncata (X95935). Among the Embryophyta reads, 2.3% 
and 40.4% in Evl2Euk and Evl3Euk respectively are highly similar 
(98-99%) to CratyUa sp. (JXl 58808) and chickpea Cicer arietinum 
(AHIIO 11 38308) that belong to Fabaceae commonly known as 
legume, the 3rd-largest land plant family and economically import- 
ant"'''. Dust events that transport pollen long-distances introduce 
vegetation changes and are prone to errors in studies (paleoecology) 
that interpret past local vegetation based on presence of pollen". 

Sequences affiliated with class Pinopsida (mostly conifers) were 
the 2nd largest group (25.4%) of the Embryophyta in Evl3Euk, 
none of them were found in Evl2Euk. This group included reads 
highly similar (99-100%) to e.g., Cupressus gigantea (EF053166), 
TetracUnis articulata (EU161293), Juniperus morrisonicola (EF67 
3744), Chamaecyparis pisifera (EF053165), Picea morrisonicola 
(AB026939), Pinus luchuensis (D38246). Finally, flowers and grass 
also appeared only in Evl 3Euk as 2.8% of the Embryophyta reads that 
are similar (97-100%) to Arabidopsis thaliana (X16077), Arctium 
lappa (JF703098), Tagetes sp. Nickrent 3061 (U42501), Sinapis alba 
(X17062), Lolium multiflorum (AY846367), Pestuca rubra 
(AF168844), and other species. Air masses originating in southern 
Europe do not usually harbor particles but include primary biological 
aerosol consisting of viruses, bacteria, fungal spores and plant poUen. 
The January dust storm here may have collected the particles in the 
Sinai desert on its way to the Negev. 

Alveolata. This superphylum is a monophyletic group of primarily 
single-celled eukaryotes that have adopted extremely diverse modes 
of nutrition such as predation, photoautotrophy and intracellular 
parasitism. Most alveolates fall into 3 main subgroups: ciliates, 
dinoflagellates and apicomplexans. Ciliates are one of the most 
important groups of protists common to water tarns and soils, and 
have many ecto- and endosymbiotic members, as well as some oblig- 
ate and opportunistic parasites. Both dust samples (Evl2Euk and 



Evl 3Euk) encompass 1.7% and 0. 1 % respectively of such reads; most 
reads from the former are identical to Halteria spp. with typical size 
of 20-50 |im and predominantly found in fresh water habitats. 

Arthropoda. Members of this phylum are exoskeleton-containing 
invertebrates with a segmented body and jointed appendages, and 
include insects, arachnids, crustaceans and myriapods. Evl2Euk 
contained 2.6% reads, most of which are affiliated with unclassified 
Metazoa/Animalia (with similarity of £93% in the BLAST database) 
and some belong to infraclass Neoptera — a group that includes most 
of the winged insects. Evl3Euk included only 0.26% such sequences. 

Unclassified eukarya. According to SILVA database classification, 
reads from Evl2Euk and Evl3Euk include respectively 47.5% and 
1.8% unclassified eukaryotes; they were therefore further analyzed 
using the BlastN program: about half (23%) of the Evl2Euk reads 
were identical to Rhizophlyctis rosea (NG017175), a soil cellulose- 
decomposing zoospores-producing fungus^'. Additional 4.5% were 
identical to other ubiquitous zoospore-producing soil-dwelling fun- 
gal species belonging to Gaertneriomyces, Spizellomyces, Powello- 
myces and other genera (e.g., JN020240, GU568157, JN940943, 
FJ827646, AY349038, FJ827660, HQ901755, HQ901759). Addi- 
tional 7.3% of Evl2Euk unclassified reads were identical to sequences 
of nematodes such as Aphelenchus avenae (AY284640 and 
AB731165), which is mycophagous and capable of withstanding 
droughts^". No nematodes were observed in Evl 3Euk. The remaining 
(12.7%) of unclassified eukaryotic reads from Evl2Euk appeared as 
singletons with low simOarity to known sequences. Unclassified 
Evl3Euk reads include only 69 sequences, 9 of which showed high 
similarity (97%) to an uncultured ciliate clone QD09 (HQ909037). 

Concluding remarks. The results demonstrate that the diversity of 
organisms in airborne dust is higher than previously reported and 
may rival that of the other (terrestrial or aquatic) environments." '^ 
Dust of Evl3, of south European origin, was rich with reads affiliated 
to Betaproteobacteria, and to Embryophyta (land plants) in general, 
particularly conifers. On the other hand, dust of Evl 2 (of north 
African origin) contained significantly more Actinobacteria, fungi, 
unclassified bacteria and eukaryotic sequences. Reads affiliated to 
allergenic and pathogenic species existed in both airborne dust 
samples. Intercontinental aerobiology studies will aid in developing 
predictive models for disease dispersaP'. Despite its limitations'", 
culturing data may be valuable by knowing what species remain 
viable after long distance atmospheric transport™. 

This study will likely affect the methodology to analyze climatic 
factors, soil sources, levels of particulate matter and air biology assoc- 
iated with dust storms, and the assessment of the possible impacts on 
the ecosystems and hazards to public health. 

Methods 

Dust sample collection and physico-chemical measurements. The samples of 
outdoor dust were collected in Be'er Sheva (Israel) on a building roof within the 
campus of Ben-Gurion University of the Negev, during storm events that occurred in 
the "dust season" (October - JMay): on December 20, 2012 (Evl2) and on January 7, 
2013 (Evl3). Both events were typically strong storms in the studied area. The spatial 
extent of each event was recorded by satellite images (MODIS Level 1 and VIIRS Level 
1) (Fig. 1). The air mass transport in the region was detected through Backward 
Trajectories model (NOAA/ARL HYSPLIT-4)'' prior and during the events. The 
wind directions at 10 m level were retrieved by the POSEIDON System. The 
following major pollutants and meteorological variables were recorded: NO2 (ppb), 
NO (ppb), CO (ppm), SO2 (ppb), O3 (ppb), wind speed (m s"'), air temperature ("C), 
and relative humidity (%). Atmospheric concentrations ((ig m"'') of PMio and PM2.5 
and total settled dust (g m"^) were recorded simultaneously during the storms in a 
monitoring system nearby the dust samplers. 

The sterilized settling-dust collectors consist of a rectangular plastic tray (40 X 25 
X 10 cm) filled with layers of glass/quartz marbles (10 mm in diameter). 
Atmospheric dust particles that cross the tray aperture are trapped in the marble 
matrix due to the matrix cohesion and roughness along with reduced air velocity near 
the surface of the tray. At the end of each event the dust samples were moved into 
sterilized glass vials for size distribution and elemental analyses. Size distributions of 
particles over the range of 0.08 to 2000 [im were obtained by a high-resolution laser 
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diffractometer (ANALYSETTE 22 MicroTec Plus). Each sample was dispersed by 
sonication {at 38 kHz) in a Na-hexametapliosphate solution (0.5%), then transferred 
to a fluid module of the instrument (containing deionized water), and subjected to 3 
consecutive 1 min runs at a medium pump speed of 6 L min"\ The data were 
processed using the Mie scattering model (RI — 1.56, AC — 0.01) with an error < 
5.0%. The MasControl software was employed to determine the following relevant 
parameters: mean size, median, modes in multiple modal distributions, sorting 
values, size fraction weights. 

Elemental analyses were performed by X-Ray Fluorescence (XRE) method using 
XRF Spectrometer, Panalytical Co., model Axios [WDXRF (wavelength dispersive), 
1 kW]. Omnian software was used for quantitative analysis. 

DNA extraction, PGR amplitication and pyrosequencing. Total DNA was 
extracted from 0.5 g of the 2 dust samples (Evl2 and Evl3) using a Power Soil 
Isolation Kit (MoBio Laboratories, Carlsbad, CA) with the bead-beating protocol 
supplied by the manufacturer, and quality and concentration were assessed using a 
NanoDrop ND 1000 spectrophotometer (NanoDrop Technologies, Wilmington, 
DE). 

Low concentration (~5 ng |jl~^) of the total genomic DNA of each sample was 
amplified using a Biometra TGradient thermocycler (Biometra, Gottingen, Germany) 
and 2 pairs of primers: the universal bacterial primers 8F, AGAGTTTGATYM- 
TGGCTCAG and 907R, CCGTCAATTCMTTTGAGTTT to amplify a -900 bp 
region of the 16S rRNA genes, and the universal eukaryotic primers NSF370/18, 
AGGGYTCGAYYCCGGAGA and NSR1787/18, CYGCAGGTTCACCTACRG tar- 
geting 18S rRNA genes amplifying the —1400 bp. The reaction mixture included 
12.5 pi ReddyMlx (ABgene, Surrey, UK), 1 pi of 10 |j,M concentrations of each 
primer (forward and reverse), 1 \x\ of 25 mM bovine serum albumin, 1-3 pi of the 
sample genomic DNA (10-40 ng i^l— 1), and water for a total volume to 25 pi. An 
initial denaturation hot start of 4 min at 95^C was followed by 30 cycles of the 
following incubation pattern: 94''C for 30 s, 54'^C for 30 s, and 72'^C for 1-1.5 min. 

The amplicons were submitted to Molecular Research laboratory (MR DNA, 
Shallowater, TX) for PCR optimization and pyrosequencing analyses utilized Roche 
454 (454 Life Sciences, Branford, CT). The company's primer 27Fmod 
AGRGTTTGATCMTGGCTCAG was used to obtain bacterial reads of the hyper- 
variable regions V1-V2 from the dust samples (Evl2 and Evl3), designated Evl2Bac 
and Evl3Bac, and the primer Euk528F CCGCGGTAATTCCAGCTC was used for 
the V4 eukaryotic reads designated Evl2Euk and Evl3Euk. 

Analysis of organismic diversity. In total, 15,019 bacterial and 15,424 eukaryotic 
sequence reads were recovered and subjected to MOTHUR software version 1.9.1^", 
leaving 4020 bacterial reads for each sample, Evl2Bac and Evl3Bac, and 3779 
eukaryotic reads for each Evl2Euk and Evl3Euk, for further considerations. 

The sequences were trimmed using MOTHUR software for removal of primers, 
barcodes, ambiguous nucleotides, long homo-polymers, reads below a minimum 
quality score of 25 and sequences shorter than 150 bases. Sequences were aligned, 
checked for chimeras, filtered, and classified using the Ribosomal Database Project II 
website (RDP Release 10) for bacteria (Table S2) and SILVA reference files for 
eukarya (Tables S3 and S4). Reads classified as Chloroplast, Mitochondria or 
"unknown" were also removed. The sub. sample function in MOTHUR, used to 
randomly select equaUy-sized reads from each library, further reduced data sets for 
subsequent analyses. The pyrosequencing population data were also analyzed by 
multiple sequence alignments using the dist.seqs function™. 

MOTHUR operational taxonomic units (OTUs) analyses were conducted (by a 
0.03 distance level) to present a-diversity: Chaol richness estimators, the inverse 
Simpson diversity index and rarefaction curves at 0.03 distances, which were plotted 
on a line chart using Microsoft Excel. PC-ORD 6 software (MJM Software Design, 
OR) was used to calculate beta- diversity. 

Raw sequencing data was deposited in the MG-RAST (metagenomics.anl.gov) 
archive under accession numbers: 4534886.3 for Evl2Bac; 4534887.3 for Evl2Euk; 
4534888.3 for Evl3Bac; 4534889.3 for Evl3Euk. 
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